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We define a natural ensemble of trace preserving, completely positive quantum maps and present 
algorithms to generate them at random. Spectral properties of the superoperator $ associated with 
a given quantum map are investigated and a quantum analogue of the Frobenius-Perron theorem 
is proved. We derive a general formula for the density of eigenvalues of $ and show the connection 
with the Ginibre ensemble of real non-symmetric random matrices. Numerical investigations of the 
spectral gap imply that a generic state of the system iterated several times by a fixed generic map 
converges exponentially to an invariant state. 
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■ Random dynamical systems are subject of a considerable scientific interest. Certain statistical properties of classical 
| deterministic systems which exhibit chaotic behavior can be described by random processes. A similar link exists also 
in the quantum theory, since quantum analogues of classically chaotic dynamical systems display spectral properties 
characteristic of ensembles of random matrices [l[ ■ 

In the case of an isolated quantum system the dynamics is represented by a unitary evolution operator which acts on 
the space of pure quantum states. To characterize such evolution operators for periodically driven, quantum chaotic 
systems the Dyson circular ensembles of random unitary matrices are often used. The symmetry properties of the 
£h ' system determine which universality class should be used 0]. 

In a more general case of a quantum system coupled with an environment one needs to work with density operators. 
In such a case any discrete dynamics is governed by a quantum operation which maps the set of all density operators 
of size N into itself. Investigations of the set Sj^f of all quantum operations play the key role in the field of quantum 
, information processing Q, since any physical transformation of a state carrying quantum information has to be 
£> ' described by an element of this set. 

To process quantum information one needs to transform a given quantum state in a controlled way. When designing 
a sequence of quantum operations, which constitutes a quantum algorithm, it is important to understand the properties 
£T) | of each operation and to estimate the influence of any imperfection in realization of an operation on the final result. 

In order to describe an effect of external noise acting on a quantum system one often uses certain models of random 
quantum operations. On the other hand, random operations are purposely applied to obtain pseudo-random quantum 
circuits @,@]. These possible applications provide motivation for research on random quantum operations. 

The main aim of this work is to construct a general class of random quantum operations and to analyze their 
properties. More formally, we introduce a natural probability measure which covers the entire set Sn of quantum 
operations and present an efficient algorithm to generate them randomly. We investigate spectral properties of 
supcroperators associated with quantum operations and infer conclusions about the convergence of any initial state 
subjected to the repeated action of a given random operation to its invariant state. It is worth to emphasize that the 
spectral properties of superoperators are already accessible experimentally, as demonstrated by Weinstein et al. Q in 
a study of a three-qubit nuclear magnetic resonance quantum information processor. 

We start reviewing the classical problem, in which classical information is encoded into a probability vector p of 
length N. The set of all classical states of size N forms the probability simplex Ajv-i- A discrete dynamics in this 
set is given by a transformation p^ = SijPj, where S is a real square matrix of size N, which is stochastic, i.e.: 

a) Sij > for i, j = 1,...,N; 

b) Eti% = 1 for all .7 = 1,..., N. 
If a real matrix S satisfies assumptions (a) and (b) then 

i) the spectrum {z.i}fL 1 of S belongs to the unit disk, \zi\ < 1, and the leading eigenvalue equals unity, z\ = 1; 

ii) the eigenspace associated with z\ contains (at least) a real eigenstate p lm , which describes the invariant measure 
ofS. 

If additionally 

c) Ef=iSij = 1 for alii = 1,...,N 
then the matrix S is called bistochastic (doubly stochastic) and 

Hi) the maximally mixed state is invariant, p>i nv = (1/N,.. ., 1/iV). 
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This is a form of the well-known Frobenius Perron (FP) theorem and its proof may be found e.g. in [1, Q. 

To generate a random stochastic matrix it is convenient to start with a square matrix X from the complex Ginibre 
ensemble, all elements of which are independent complex random Gaussian variables. Then the random matrix 

N 

Sa ■■= 1^-17 E \ x ^\ 2 > (!) 

m— 1 

is stochastic [Io| , and each of its columns forms an independent random vector distributed uniformly in the probability 
simplex Ajv-i- 

Let us now discuss the quantum case, in which any state can be described by means of a positive, normalized 
operator, p > 0, Trp = 1. Let M.n denote the set of all normalized states which act on a iV-dimensional Hilbert space 
Hjy. A quantum analogue of a stochastic matrix is given by a linear quantum map $ : Aijy — > A4n which preserves 
the trace, and is completely positive (i.e. the extended map, <&<g> Ik, is positive for any size K of the extension Ik)- 
Such a transformation is called quantum operation or stochastic map, and can be described by a matrix $ of size iV 2 , 

p' = $/9 Or p^u = Qm^Pnv ■ (2) 

nv 

It is convenient to reorder elements of this matrix pH defining the so-called dynamical matrix D, 

£>($) EE $ fl SO that D mn = $m» , (3) 

nv nv 

since the matrix D is Hermitian if the map <E> preserves Hermiticity. Note that the above equation can be interpreted 
as a definition of the operation of reshuffling (realignment) of a matrix X, denoted by X R . This definition is 
representation dependent, and corresponds to exchanging the position of some elements of a matrix. Obviously 
reshuffling is an involution: D R = $. 

Due to the theorem of Choi [l2T |. a map $ is completely positive (CP) if and only if the corresponding dynamical 
matrix D is positive, D > 0. Hence D can be interpreted as a positive operator actin g o n a composed Hilbert space 
H := TLa 8%. Any completely positive map can also be written in the Kraus form [13j, 

p^ p' = Y,A iP A\. (4) 

i 

The set of Kraus operators A; allows one to write down the linear superoperator as 

$ = E A i ® Ai , (5) 

i 

where Ai is the complex conjugate of Ai, while A\ denotes the adjoint operator: A\ = Af . 

The map $ is trace preserving, Tip' = Tip = 1, if ^2iA\Ai = Ijv. This condition is equivalent to a partial trace 
condition imposed on the dynamical matrix, 

Tr A D($) = t N , (6) 

which implies Tr D = N. 

Since the dynamical map of an operation $ is positive and normalized, the rcscalcd matrix D/N may be considered 
as a state in an extended Hilbert space Ha ® of size N 2 . Stochastic maps and states on the extended space 
Ha ®Hb are related by the so-called Jamiolkowski isomorphism 0, [l4| • Making use of the maximally entangled 
bipartite state \ip + ) = ^7=Xa=i |M) > tms can expressed as D($)/N = ($ (g> 1)|V'+)(V'+I ■ For completeness we 
sketch here a compact proof of the latter fact. 

Proof. The projector \il> + )(ip + \ can be recastcd in the form X)fj=i N)0U ® l*)0'|s , in which the matrix elements 
on Ha and Hb are explicit. Then, (§a <8> 1b)|V'+)(V , + I = ^ j=i ^ (K)01)y4 ® (N)01)b > wn ich in turn can be 
expressed in coordinates as -4 X^f^i H i/=i < ^ > ' Jly ^^(^jl ■ Thus the result follows from definition . ■ 

' ' ' ij 

A quantum map is called unital if it leaves the maximally mixed state invariant. It is so if the Kraus operators 
satisfy J2iAiA\ = Ijv- The unitality condition may also be written in a form Tr^-D = Ijv, dual to . A CP 
quantum map which is trace preserving and unital is called bistochastic. 

Spectral properties of positive operators were studied in the mathematical literature [l5l,[l(| for a general framework 
of C*-algebras. Here we analyze spectral properties of the operator $ corresponding to a stochastic map and formulate 
a quantum analogue of the Frobenius Perron theorem. 
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Let $ be a complex square matrix of size N 2 , so that it represents an operator acting in a composite Hilbcrt space 
TL N 2 ~ Ha ®Ti.B- Let us order its complex eigenvalues according to their moduli, \zi\ > \z 2 \ > ■ ■ ■ > \z N 2\ > 0. 
Assume that <I> represents a stochastic quantum map, i.e. it satisfies 
a') <P R > 0; b') Y, k * fc fc = % so that Tr A $ fl = 1, 

»i 

where the reshuffling operation, denoted by R , is dehned in eq. ([3]) . Then 

i') the spectrum {zi]f =l of $ belongs to the unit disk, \zi\ < 1, and the leading eigenvalue equals unity, z\ = 1, 

if) one of the corresponding eigenstates forms a matrix u> of size N which is positive, normalized (Tr ui = 1) and is 

invariant under the action of the map, = to. 

If additionally 

d) Ty b <§> r ; = 1, then the map is called bistochastic and 
Hi') the maximally mixed state is invariant, to = 1/JV. 

Proof. Assumption (a') implies that the quantum map is completely positive while (b') implies the trace preserving 
property. Hence $ is a linear ma p w hich sends the convex compact set M. N of mixed states into itself. Due to 
the Schauder fixed-point theorem [17| such a transformation has a fixed point — an invariant state uj > such that 
$w = oj. Thus z\ = 1 and all eigenvalues fulfil | < 1, since otherwise the assumption that $ maps the compact set 
Mn into itself would be violated. □ 
The spectral properties of stochastic maps are similar to those of classical stochastic matrices discussed in [l8[ . As 
noted in [lj| Ho] the spectrum of the superoperator is symmetric with respect to the real axis - see Fig . 1. Such 
spectra for a quantum map corresponding to a classically chaotic irreversible system where studied in [211 ] . 




FIG. 1: Spectra of an ensemble of random superoperators a) 800 operations for N = 2 and b) 300 operations for iV = 3. 

Let us analyze briefly some particular cases of quantum operations. For any unitary rotation. $(p) = UpW the 
Kraus form ^ consists of a single operator only A\ = U . Thus according to ^ the superoperator is given by a 
unitary matrix $ = U <E> U of size TV 2 . Denote the eigenphases of U by cti where i = 1, . . . N. Then the spectrum of 
<i> consists of TV 2 phases given by at — ctj for i, j = 1, . . . ,N. All diagonal phases for i = j are equal to zero, hence 
the leading eigenvalue of the superoperator $ , Z\ = 1 , exhibits multiplicity not smaller than N. 

Consider now a quantum map for which $> R is diagonal. This special case can be treated as classical, since then 
$ describes a classical dynamics in Ajy-i, while the generalized quantum version of the FP theorem reduces to its 
standard version. Reshaping a diagonal dynamical matrix of size ./V 2 one obtains then a matrix S of size N, where 
Sij = (no summation performed!). Then assumption (a') (all diagonal elements of $> R are positive) gives (a), 
ij 

while the trace preserving condition (b') implies the probability preserving condition (b). Similarly, the additional 
condition (c') for quantum unitality gives condition (c) which imposes that the uniform vector is invariant under 
multiplication by a bistochastic matrix S. Similarly, conclusions (i'), (ii') and (iii') of the quantum version of the 
theorem imply conclusions (i-iii) of the standard (classical) Frobcnius-Pcrron theorem. 

We are interested in defining an ensemble of random operations Q • A simple choice of random external fields [22l | , 
defined as a convex combination of an arbitrary number k of unitary transformations, p' = 'Yli = iPiVipV i , produces 
bistochastic operations only. 

Let us then consider first a method of constructing random states by N x M rectangular random complex matrices 
of the Ginibre ensemble |23 ; ]. Taking 



p := XX^/TrXX^ , 



(7) 
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we get a positive normalized state p. For M = 1 we obtain a recipe to generate random pure states, while for M = N 
the measure induced by the Ginibre ensemble coincides with the Hilbert-Schmidt measure in A4n [24[ and the average 
purity, (Trp 2 ) H s, scales as 1/N 

Here we propose an analogous algorithm of constructing a random operation: 

1) fix M > 1 and take a N 2 x M random complex Ginibre matrix X; 

2) find the positive matrix Y := TtaXX^ and its square root y/Y ; 

3) write the dynamical matrix (Choi matrix) 




4) reshuffle the Choi matrix according to ([3]) to obtain the superoperator <E> = D R , and use it as in @ to produce a 
random map. 

It is not difficult to check that the relation © holds due to ([5]) , so the random map preserves the trace. Such 
a renormalization to obtain the Choi matrix was independently used in (25[. This method is simple to apply for 
numerical simulations, and exemplary spectra obtained in the case M = N 2 are shown in Fig.l. 

For larger N, the subleading eigenvalue modulus r = \z2\ is smaller, so the convergence rate of any initial po to 
the invariant state uj occurs faster. To demonstrate this effect we studied the decrease of an average trace distance 
in time, L(t) = (Tr|$ (pa) — where the average is performed over an ensemble of initially pure random states, 

Pa = IV 1 ) ("01- Numerical results confirm an exponential convergence, L(t) ~ exp(— at). The mean convergence rate 
(a}$, averaged over an ensemble of random operations, increases with the dimension N like logiV, with slope very 
close to unity - see Fig. 2. 




FIG. 2: a) Average trace distance of random pure states to the invariant state of $ as a function of time for 
N = 4(»), 6(B), 8(*); b) mean convergence rate (a)* as a function of the system size N, plotted in log scale. 



To explain these findings we need to analyze spectral properties of an ensemble of random operations. Consider 
first a random Choi matrix D obtained by the above algorithm from a Ginibre matrix X, generated according to the 
distribution oc exp(— TiXX^). Then D is of the Wishart type and has the distribution 

P(D) oc J dX cxp(-TrX X r )5(D - WXX^W^) , (9) 

where W := t <g> (TtaXX^)- 1 / 2 . This integral can be rewritten with help of another Delta function, 

P(D) oc fdYjdX exp(-TrY)5(D - WXX t W t )S(Y -Tr A XX^) . (10) 

Using the Delta function property and taking the Jacobians into account we arrive at 

P(D) oc Act(D M - N2 ) 8(Tr A D -1), (II) 

which shows that there are no other constraints on the distribution of the Choi matrix, besides the partial trace 
condition and positivity. This is equivalent to saying that the M Kraus matrices Ai, which form a map <&, constitute 
a (MN) x N truncated part of a unitary matrix U of size NM . A natural assumption is that the matrix U is 
distributed according to the Haar measure. Unitarity constraints become weak for large N, so the non-Hermitian 
N x N truncations Ai are described by the complex Ginibre ensemble [23| : Their spectra cover uniformly the disk of 
radius 1/y/N in the complex plane [2r| . 
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Therefore we are in position to present an alternative algorithm of generating the same ensemble of random oper- 
ations, which has a simple physical interpretation: 

(1') Choose a random unitary matrix U according to the Haar measure on U(NM); 
(2') Construct a random map defined by 

p' = Tt M [U{p® \v)(v\)rf\ . (12) 

Hence this random operation corresponds to an interaction with an A/-dimensional environment, initially in a random 
pure state \ v). Of a special importance is the case M = N 2 , for which the term with the determinant in (jTTJ) disappears, 
so the matrices D are generated according to the measure analogous to the Hilbcrt-Schmidt measure. This case 
provides thus generic dynamical matrices of a full rank and can be recommended for numerical implementation. 

To analyze spectral properties of a supcropcrator $ let us use the Bloch representation of a state p, 

N 2 -l 



i=0 



where A* are generators of SU(N) such that tr(A l A-?) = <5 y and A = l/y/N. Since p — p* , the generalized Bloch 
vector t = [to, . . . , Tjv^-i], also called coherence vector, is real. Thus the action of the map $ can be represented as 

t' = $(r) = Ct + k, (14) 

where G is a real asymmetric contraction matrix of size iV 2 — 1 while k is a translation vector, which vanishes for 
bistochastic maps. Their elements can be expressed in terms of the Kraus operators, e.g. Cy = Tr^ fe AMfeA^Aj., 
while Ki = TrA'A TkA^. Thus there exists a real representation of the superoperator 



1 

K C 



(15) 



Eigenvalues of G, denoted by {A*}^ -1 , are also eigenvalues of We are going to study the case M = N 2 , for which 
the distribution (jTTJ) simplifies. Like for the real Ginibre ensemble [13, [H[ one may derive in this case the measure in 
the space of eigenvalues Aj of G 

dp(A) = \Y[dAi\ JJ |A fe - Ai\ G(A) , (16) 

i k<l 

where G(A) is given by the distribution of (real) traces, 

N 2 -l N 2 -l 



u=l i=l 



G(A) :=([] d^W-TrCr)) (17) 

with the average 

(f(C)) := I dCdn Q(D > 0) f(C) . 



The domain of integration is given by the conditions for complete positivity, D > 0, which is not easy to work with, 
even for N = 2. For large N we can expect that these conditions do not play an important role, so the dependence 
G(A) is weak, and the measure for G can be described by the real Ginibre ensemble of non-symmetric Gaussian 
matrices. The spectrum of such random matrices consists of a component on the real axis, the probability density 
of which is given asymptotically by the step function P(x) = ^0(1 — |x|) (29l. [30|. while remaining eigenvalues cover 
uniformly the unit circle according to the Girko distribution 

To analyze the spectra of random operators $ one needs to set the scale. The mean purity of a random state a 
of size A^ 2 behaves as N~ 2 [l(| and D = Na, thus the average Tr£> 2 = Tr$$^ is of the order of unity. Hence, the 
rescaled matrix <£>' := A^$ of size K = N 2 has the normalization Tr$'($')^ sw K, which assures that the radius of the 
circle is equal to unity. 

Thus we arrive at the following conjecture: for large N the statistical properties of a rescaled random superoperator 
N<f> are described by the real Ginibre ensemble. We confirmed this conjecture by a detailed numerical investigation. 
Fig. 3 shows the density of complex eigenvalues of random superoperators for A^ = 10 with AI = N 2 and the 



FIG. 3: a) Distribution of complex eigenvalues of 10 4 rescaled random operators $' already for TV = 10 can be approximated 
by the circle law. b) Distribution of real eigenvalues P(x) of $' plotted for TV = 2, 3, 7 and 14 tends to the step function, 
characteristic of real Ginibre ensemble. 



distribution P(x) of the real eigenvalues. As the spectrum of the rescaled operator $' = TV$ tends to be localized in 
the unit circle, we infer that the size of the sublcading eigenvalue r = \z2\ of $ behaves as 1/TV, hence the convergence 
rate a scales as In TV. 

Numerical studies were also performed for random maps acting on states of a fixed dimension TV. In this case the 
subleading eigenvalue of a random map (|12p decreases with the varying size of the environment M as r ~ 1/ vM '. 
Similar investigations were also performed under a constraint that the dynamical matrix D is diagonal. In this case the 
assumption M = 1 allows to obtain a random stochastic matrix S of size TV, such that each of its columns is generated 
independently with respect to the flat measure in the (TV — 1) dimensional simplex of probability distributions. 
Analyzing the average trace of SS T we infer that in this case the complex spectrum can be described by the Girko 
distribution, which covers uniformly the disk of radius r ~ l/\/~N. These spectral properties of random stochastic 
matrices, confirmed by our numerical results, were rigorously analyzed in a recent paper of Horvat pjl] ]. 

In this work we analyzed superoperators associated with quantum stochastic maps and their spectral properties 
and formulated a quantum analogue of the Frobenius-Perron theorem. We defined an ensemble of random operations, 
presented an explicit algorithm to generate them, and showed an exponential convergence of a generic state of the 
system to the invariant state under subsequent action of a fixed map. We demonstrated that for a large dimension 
of the Hilbcrt space, used to describe quantum dynamics, the spectral properties of a generic superoperator can be 
described by the Ginibre ensemble of real random matrices. 

It is a pleasure to thank M.D. Choi, W. Slomczyhski and J. Zakrzcwski for helpful remarks. We acknowledge 
financial support by the SFB Transregio-12 project, the European grant COCOS, and a grant number N202 099 
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